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This article reviews and updates the Standard Model prediction of the r lepton g— 2. 
Updated QED and electroweak contributions are presented, together with new values of 
the leading-order hadronic term, based on the recent low energy e+e — data from BaBar, 
CMD-2, KLOE and SND, and of the hadronic light-by-light contribution. The total 
prediction is confronted to the available experimental bounds on the r lepton anomaly, 
and prospects for its future measurements are briefly discussed. 

1. Introduction 

Numerous precision tests of the Standard Model (SM) and searches for its possible 
violation have been performed in the last few decades, serving as an invaluable 
tool to test the theory at the quantum level. They have also provided stringent 
constraints on many "New Physics" (NP) scenarios. A typical example is given by 
the measurements of the anomalous magnetic moment of the electron and the muon, 
where recent experiments reached the fabulous relative precision of 0.7 ppb^and 0.5 
ppmpl respectively. These experiments measure the so-called gyromagnetic factor 
<7, defined by the relation between the particle's spin s and its magnetic moment /2, 

where e and m are the charge and mass of the particle. In the Dirac theory of 
a charged point-like spin-1/2 particle, g = 2. Quantum Electrodynamics (QED) 
predicts deviations from Dirac's value, as the charged particle can emit and reabsorb 
virtual photons. These QED effects slightly increase the g value. It is conventional 
to express the difference of g from 2 in terms of the value of the so-called anomalous 
magnetic moment, a dimensionless quantity defined as a = {g — 2)/2. 

The anomalous magnetic moment of the electron, a e , is rather insensitive to 
strong and weak interactions, hence providing a stringent test of QED and lead ing 
to the most precise determination of the fine-structure constant a to datePEl On 
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the other hand, the g—2 of the muon, a M , allows to test the entire SM, as each of 
its sectors contributes in a significant way to the total prediction. Compared with 
a e , a M is also much better suited to unveil or constrain NP effects. Indeed, for a 
lepton I, their contribution to a; is generally expected to be proportional to mf /A 2 , 
where m; is the mass of the lepton and A is the scale of NP, thus leading to an 
(m M /m e ) 2 ~4x 10 4 relative enhancement of the sensitivity of the muon versus the 
electron anomalous magnetic moment. This more than compensates the much higher 
accuracy with which the g factor of the latter is known. The anomalous magnetic 
moment of the r lepton, a T , would suit even better; however, its direct experimental 
measurement is prevented by the relatively short lifetime of this lepton, at least at 
present. The existing limits are based on the precise measurements of the total and 
differential cross sections of the reactions e + e~ — > e + e - T + T~ and e + e _ — > Z — > 
t + t~7 at LEP energies. The most stringent limit, —0.052 < a T < 0.013 at 95% 
confidence level, was set by the DELPHI collaboration,^ and is still more than an 
order of magnitude worse than that required to determine a T . 

In the 1990s it became clear that the accuracy of the theoretical prediction of 
the muon g—2, challenged by the E821 experiment underway at Brookhavenpl was 
going to be restricted by our knowledge of its hadronic contribution. This problem 
has been solved by the impressive experiments at low-energy e + e~ colliders, where 
the total hadronic cross section (as well as exclusive ones) were measured with high 
precision, allowing a significant improvement of the uncertainty of the leading-order 
hadronic contribution!^^ As a result, the accuracy of the SM prediction for et M now 
matches that of its measurement. In parallel to these efforts, very many improve- 
ments of all other sectors of the SM prediction were carried on by a large number of 
theorists (see Refs. [8J H2 for reviews). All these experimental and theoretical devel- 
opments allow to significantly improve the theoretical prediction for the anomalous 
magnetic moment of t lepton as well. 

In this article we review and update the SM prediction of a T , analyzing in detail 
the three contributions into which it is usually split: QED, electroweak (EW) and 
hadronic. Updated QED and EW contributions are presented in Sees. [2] and [3J 
new values of the leading-order hadronic term, based on the recent low energy 
e+e~ data from BaBar, CMD-2, KLOE and SND, and of the hadronic light-by-light 
contribution are presented in Sec. [4] The total SM prediction is confronted to the 
available experimental bounds on the r lepton g—2 in Sec. [5) and prospects for its 
future measurements are briefly discussed in Sec. [H where conclusions are drawn. 



2. QED Contribution to a T 

The QED part of the anomalous magnetic moment of the r lepton arises from the 
subset of SM diagrams containing only lep tons and photons. This dimensionless 
quantity can be cast in the general formM^ 

<fT =A 1+ aJ^) + aJ^)+A 3 (^^), (2) 
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where m e , and m T are the electron, muon and r lepton masses, respectively. The 
term A\, arising from diagrams containing only photons and r leptons, is mass and 
flavor independent. In contrast, the terms Ai and A3 are functions of the indicated 
mass ratios, and are generated by graphs containing also electrons and/or muons. 
The functions Aj (i = 1, 2, 3) can be expanded as power series in a/n and computed 
order-by-order 

(£ )+ , f(3 %4.) ( £)» + , f (£)%.... (3) 

Only one diagram is involved in the evaluation of the lowest-order (first-order in 

a, second-order in the electric charge) contribution - it provides the famous result 

by Schwinger A^ } = l/2P]The mass-dependent coefficients A2 and A3, discussed 

-1 01 

below, are of higher order. They were derived using the latest CODATA-^ recom- 
mended mass ratios: 

m T /m e = 3477.48(57) (4) 
m T /m AI = 16.8183(27). (5) 

The value for m T adopted by CODATA in Ref. QH m T = 1776.99 (29) MeV, is 
based on the PDG 2002 result ff^ It remained unchanged until very recently (see 
Refs. [HI H5l). when preliminary results of two new measurements (from the BellJ^ 
and KEDFiP^ detectors) were reported. The central values of the new mass values 
are slightly lower than the current world average value, but agree with it within the 
uncertainties, which are approaching that of the world average value (used in this 
work). 



2.1. Two-loop Corrections 

Seven diagrams contribute to the fourth-order coefficient A^\ one to A^' ^ {m T / 'm e ) 
and one to 4 4 '( m r/ m /i)- They are depicted in Fig. [T] As there are no two- loop 
diagrams contributing to a!j? ED that contain both virtual electrons and muons, 
A^\m T /m e ,m T /m ( ,) = 0. The mass-independent coefficient has been known for 
almost fifty years 

„ M1 197 7T 2 3_ 7T 2 , „ 

AT' = 1 h -C(3 n2 

1 144 12 4 SW 2 

= -0.328 478 965 579193 78..., (6) 

where £(s) is the Riemann zeta function of argument s. For I = e, fj, or r, the 
coefficient of the two-loop mass-dependent contribution to a^ ED , A% (1/x), with 
x = rrij/mi, is generated by the diagram with a vacuum polarization subgraph 
containing the virtual lepton j. This coefficient was first computed in the late 1950s 
for the muon g— 2 with x = m e /m M <C 1, neglecting terms of 0(x)P^ The exact 
expression for < x < 1 was reported by Elend in 1966p^ However, its numerical 
evaluation was considered tricky because of large cancellations and difficulties in the 
estimate of the accuracy of the results, so that common practice was to use series 
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Fig. 1. The QED diagrams contributing to the r lepton g— 2 at order a 2 . The mirror reflections 
(not shown) of the third and fourth diagrams must be included as well. 



expansions instead t HIHEll Taking advantage of the properties of the dilogarithm 
Li2(z) — — J* dtln(l — t) I tj^l the exact result was cast in Ref. [Din a very simple 
and compact analytic form, valid, contrary to the one in Ref. [501 a l so f° r x > 1 (the 
case relevant to a^ ED and part of aJJ BD ): 

25 In a; 



4 4) (l/x) 



36 



+ x 2 (4 + 3 lux) + - (l - 5x 2 ) x 
Li 2 (x) + Li 2 (— x) 
-U 2 {x 2 ) 







r l-x\ 


X 


In x In 1 






2 \ 


y l + x) 



2 In x In I — 

3 \ a; 



(7) 



For x = 1, Eq. |[7J) gives (1) = 119/36 — 7r 2 /3; of course, this contribution is 
already part of in Eq. |6]). Numerical evaluation of Eq. (O with the mass ratios 
given in Eqs. (HJ)-© yields the two-loop mass-dependent QED contributions to the 
anomalous magnetic moment of the t leptonp^ 

A ( 2 4 \m T /m e ) = 2.024284(55), (8) 
A ( 2 4) (m r /m M ) = 0.361 652 (38). (9) 

These two values are very similar to those computed via a dispersive integral in 
Ref. [25] (which, however, contain no estimates of the uncertainties). Equations ([8]) 
and © are also in agreement (but more accurate) with those of Ref. Adding up 



Eqs. 



and (J9J) one gets J 



A 



C ( / } = 2.057457(93) 



(10) 



(note that the uncertainties in m T /m e and m T /m^ are correlated). The resulting 
error 9.3 x 10^ 5 leads to a 5 x 10~ 10 uncertainty in a^ ED . 
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2.2. Three-loop Corrections 

More than one hundred diagrams are involved in the evaluation of the three-loop 
(sixth-order) QED contribution. Their analytic computation required approximately 
three decades, ending in the late 1990s. The coefficient A^ arises from 72 diagra ms. 
Its exact expression, mainly due to Remiddi and his collaborators, r 

eadsESEII 

A W M * _ 215 _ _239_ 4 28259 
1 " 72^ C(3) 24 C(5) 2160^ + 5184 + 
139 _ 298 2l „ 17101 , 



-C(3)-— 7r 2 ln2 



-7T" 



18 ,w 9 810 
100 
+ - 

1.181241456 587.... (11) 



Li 4 (l/2) + — (In 2 2 - 7T 2 ) In 2 2 



This value is in very good agreement with previous results obtained with numerical 
methods 

The calculation of the exact expression for the coefficient A^' (mi/mj) for ar- 
bitrary values of the mass ratio mi/rrij was completed in 1993 by Laporta and 
Remiddi ^^ l (earlier works include Refs.[3"Tj). Let us focus on a^ ED (I = r, j = e, fx). 
This coefficient can be further split into two parts: the first one, A| (m;/mj,vac), 
receives contributions from 36 diagrams containing either electron or muon vac- 
uum polarization loops whereas the second one, (mi/rrij, Vol), is due to 12 
light-by-light scattering diagrams with either electron or muon loops. 1 ^ 1 The ex- 
act expressions for these coefficients are rather complicated, containing hundreds 
of polylogarithmic functions up to fifth degree (for the light-by-light diagrams) and 
complex arguments (for the vacuum polarization ones) - they also involve harmonic 
polylogarithmsP2l Series expansions were provided in Ref. [30] for the cases of phys- 
ical relevance. 

Using the recommended mass ratios given in Eqs. ([4} and ([5]), the following 
values were recently computed from the full analytic expressions!^ 

Af ] (m T /m e , vac) = 7.256 99 (41) (12) 

A { 2 6) (m T /m e , lbl) = 39.1351 (11) (13) 

Af ] (m T /m M , vac) = -0.023 554 (51) (14) 

A 2 6) (m T /r7i M ,lbl) = 7.033 76 (71). (15) 

Almost identical values were obtained employing the approximate series expansions 
of Ref. 30: 7.25699(41), 39.1351(11), -0.023564(51), 7.03375(71) W The previous es- 
timates of Ref. [251 were different: 10.0002, 39.5217, 2.9340, and 4.4412 (no error esti- 
mates were provided), respectively; they are superseded by the results in Eqs. (fl2|) - 
(fTS"]) . derived from the exact analytic expressions. The estimates of Ref. l33l compare 
slightly better: 7.2670, 39.6, —0.1222, 4.47 (no errors provided). In the specific case 
of A^ (m T /m^, Vol), the values of Refs . l25l and l33l differ from Eq. (fTS"]) because their 
derivations did not include terms of 0{m^/m T ), which turn out to be unexpectedly 
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large. The sums of Eqs. ([T2 ]) -([I5 j) and |[H [) -([T5 |) are 

Af\m T /m e ) = 46.3921 (15), (16) 

A^\m T /m^) = 7.010 21 (76). (17) 

The contribution of the three-loop diagrams with both electron- and muon-loop 
insertions in the photon propagator was calculated numerically from the integral 
expressions of Ref. [2TJ obtaining^! 

Af\m T /m e ,m T /m^) = 3.34797(41). (18) 

This value disagrees with the results of Refs.[33|(i-679) andl25l(2. 75316). Combining 
the three-loop results of Eqs. HI]), (16]), (HI) and OS) one finds the sixth-order QED 
coefficient J21 

Cf } = 57.9315 (27). (19) 

The error 2.7xl0 -3 induces a 3xl0 -11 uncertainty in a^ ED . The order of magnitude 
of the three-loop contribution to a^ ED , dominated by the mass-dependent terms, is 
comparable to that of EW and hadronic effects (see later). 

Contrary to the case of the elect ron an d muon g—2, QED contributions of or- 
der higher than three are not knownPSEH (An exception is the mass- and flavor- 
independent term A 4 j"^ which is however expected to be a very small part of the 
complete four-loop contribution.) Adding up all the above contributions and using 
the new value ofa derived in Refs.ElandlUa- 1 = 137.035 999 709(96), one obtains 
the total QED contribution to the g—2 of the r lepton,Sl 

a QED = 11732 4(2) x 10 -8 . (20) 

The error 8a^ ED is the uncertainty 5C? ] {a/n) 4 ~ tt 2 In 2 (m T / 'm e ){a / 'tt) 4 - 2 x 10~ 8 
assigned to a^ ED for uncalculated four-loop contributions. As we mentioned earlier, 
the errors due to the uncertainties of the 0(a 2 ) and 0(a 3 ) terms are negligible. The 
error induced by the uncertainty of a is only 8 x 10~ 13 (and thus totally negligible). 



3. Electroweak Contribution to a T 

With respect to Schwinger's contribution, the EW correction to the anomalous 
magnetic moment of the r lepton is suppressed by the ratio (m T /M w ) 2 , where M w 
is the mass of the W boson. Numerically, this contribution is of the same order of 
magnitude as the three-loop QED one. 



3.1. One-loop Contribution 



The analytic expression for the one- loop EW contribution to a Tl due to the diagrams 
in Fig. El isP3 



'(1 loop) 



5G M m 2 
24V2V 



l + -(l-4sin 2 
5 



O 



M 2 

Z.W,H 



(21) 
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where = 1.16637(1) x 10~ 5 GeV 2 is the Fermi coupling constant p3 M z , M w 
and M H are the masses of the Z, W and Higgs bosons, and 8 W is the weak mix- 
ing angle. Closed analytic expressions for a^ w (l loop) taking exactly into account 
the m^./Mg dependence (B = Z,W, Higgs, or other hypothetical bosons) can be 
found in Refs.[371 Following Ref. [351 we employ for sin 2 fly the on-shell definition}^ 
sin 2 w = 1 - Ml/Ml, where M z = 91.1876(21) GeV^and M w is the theoretical 
SM prediction of the W mass. The latter can be easily derived from the simple 
analytic formulae of Ref. [40] (see also Refs. |4T|) . 



Ma 



80.3676 - 0.05738 In 



100 GeV 



0.00892 W 



100 GeV 



GeV, (22) 



(on-shell scheme II with Aa|f } = 0.02758(35),^ a s (M z ) = 0.118(2)JI^and M top = 
171.4(2.1) GeV,S3) leading to M w = 80.343 GeV for M H = 150 GeV. This result 
should be compared with the direct experimental value M w — 80.403(29) GeVp^ 
which corresponds to a very small M H . In any case, these shifts in the M w prediction 
induced by the variation of M H from 114.4 GeV, the current lower bound at 95% 
confidence levelp^up to a few hundred GeV, only change a^ w (l loop) by amounts of 
O(10- 10 ). From Eq. ([2l]), including the tiny 0{m 2 T /M 2 z WH ) corrections of Refs.EH 
for M H = 150 GeV we get 



<^ w (l loop) 



55.1(1) x 10" 



(23) 



The uncertainty encompasses the shifts induced by variations of M H from 114 GeV 
up to a few hundred GeV, and the tiny uncertainty due to the error in m T . 







Fig. 2. One-loop electroweak contributions to a T . The diagram with a W and a Goldstone boson 
(0) must be counted twice. 



The estimate of the EW contribution in Ref.[33l a™ = 55.60(2) x 10~ 8 , obtained 
from the one- loop formula (without the small corrections of order m 2 /M 2 WH ), is 
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similar to our value in Eq. (|23|) . However, its uncertainty (2 x 10 -10 ) is too small, 
and it doesn't contain the two-loop contribution which, as we'll discuss in the next 
section, is not negligible. 

3.2. Two-loop Contribution 

The two-loop EW contributi ons to a/ (l — e,u or r) were computed in 1995 by Czar- 
necki, Krause and Marciano.™This remarkable calculation leads to a significant 
reduction of the one-loop prediction. Naively one would expect the two-loop EW 
contribution ap w (2 loop) to be of order (a/V ) x a"(l loop), but this turns out not 
to be so. As first noticed in the early 1990sp^ af w (2 loop) is actually quite sub- 
stantial because of the appearance of terms enhanced by a factor of ln(M zw /m/), 
where rrif is a fermion mass scale much smaller than M w . 

The two-loop contribution to a™ involves 1678 diagrams in the linear 't Hooft- 
Feynman gauge^ (as a check, the authors of Refs. [45] and 06] employed both this 
gauge and a nonlinear one in which the vertex of the photon, the W and the un- 
physical charged scalar vanishes). It can be divided into fermionic and bosonic 
parts; the former, a^ w (2 loop fer), includes all two- loop EW corrections containing 
closed fermion loops, whereas all other contributions are grouped into the latter, 
a^ w (2 loop bos). The expressions of Ref. [H3 for the bosonic part were obtained 
in the approximation M H 3> M W Zl computing the first two terms in the expan- 
sion in M? v z/Mg, and expanding in sin W , keeping the first four terms in this 
expansion (this number of powers is sufficient to obtain an exact coefficient of 
the large logarithms \n(M?,JmD). Recent analyses of the EW bosonic corrections 
of the g— 2 of the muoil^ relaxed these approximations, providing analytic re- 
sults valid also for a light Higgs. Considering the present M H > 114.4 GeV lower 
bound,S3! we can safely employ the results of Ref. [46] obtaining, for M H = 150 GeV 
< w (2 loop bos) = -3.06 x 10~ 8 . The neglected terms are of O(10" 9 ). 

The fermionic part of a^ w (2 loop) contains the contribution of diagrams with 
light quarks; they involve long-distance QCD for which perturbation theory can- 
not be employed. In particular, these hadronic uncertainties arise from two types 
of two-loop diagrams: the hadronic photon-Z mixing, and quark triangle loops 
with the external photon, a virtual photon and a Z attached to them (see Fig. [3]). 
The hadronic uncertainties mainly arise from the latter ones. Two approaches were 
suggested for their study: in Ref. [45] the nonperturbative effects where modeled in- 
troducing effective quark masses as a simple way to account for strong interactions. 
In view of the high experimental precision of the g— 2 of the muon, a more realistic 
treatment of the relevant hadronic dynamics was introduced in Ref. [50] within a 
low-energy effective field theory approach, later on developed in the detailed anal- 
yses of Refs. [38] [51] [52] However, from a numerical point of view, the discrepancy 
between the results provided by these two different approaches turns out to be ir- 
relevant for the present interpretation of the experimental result of the muon g— 2, 
in spite of its precision. The use of effective quark masses for the study of the g— 2 
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Hadrons 



■r 




Hadrons 



Fig. 3. Some of the fermion-loop diagrams contributing to the r anomalous magnetic moment. 



of the t - whose experimental precision is a far cry from that of the muon! - thus 
appears to be sufficient at present. The tiny hadronic 7-Z mixing terms can be eval- 
uated either in the free quark approximation or via a dispersion relation using data 
from e + e _ annihilation into hadrons; the difference was shown to be numerically 
insignificant PS! 

References [45] and |47| contain simple approximate expressions for the contribu- 
tions of the diagrams with fcrmion triangle loops shown in Fig. [3] (right). In general, 
for a lepton l — e, /1 or t, neglecting small mass-ratios, they are 



' (2 loop fer; triangle loops) 



E 



5G M mf 
24V2tt 2 



a 18 



N f llQ 2 f AC(f) 



(24) 



where l£ is the third component of the weak isospin of the fermion / in the loop, 
Qf is its charge, Nf its number of colors (3 for quarks, 1 for leptons), and 



AC(/) 



ln(M|/m 2 ) - 
ln(M|/m 2 )- 
\n(M 2 /m 2 f ) 
[(M 2 /M 2 op )/6 



5/6 

8^ 2 /27- 
2 



11/18 



(5/18)(M z 2 /M 2 op ) 



2/3]ln(Af 2 op /M 2 ) 
4/3 



if 
if 
if 

if 



m f < mi < M z 
rut = mi <C M z 
mi « m/ « M z 



(25) 



top ■ 



The contribution of the top-quark triangle loop diagram of Fig. [3] (right) with the 
Z boson replaced by the neutral Goldstone boson (0°) has also been included in 
this expression!^ It is clear from Eqs. (f24|) -([2"5" |) that the logarithms ln(A/ z ) cancel 
in sums over all fermions of a given generation, as long as to/ <C M ? , due to the no- 
anomaly condition J2f Nfl^Qj = valid within every generation!^^ This does 
not occur for the third generation due to the large mass of the top quark. Note that 
this short-distance cancellation does not get modified by strong interaction effects 
on the quark triangle diagrams (SEES 

Contrary to the case of the muon g— 2, where all fermion masses, with the excep- 
tion of m e , enter in a™ (2 loop fer; triangle loops), the approximate expressions in 
Eq. (|25|) show that this is not the case for the g—2 of the r lepton. Indeed, due to the 
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high infrared cut-off set by m T , Eq. (|24p for 1 = t does not depend on any fermion 
mass lighter than m T ; apart from m T , it only depends on M top and nib, the masses 
of the top and bottom quarks (assuming m c <m T ). The charm contribution requires 
some care, as the crude approximation provided by Eq. (|25[) for a charm lighter than 
the r lepton is valid only if m c <C m T . Clearly, this is not a good approximation, and 
the spurious shift induced by Eq. ([25)1 when m c is varied across the ro T thr eshold 
is of O(10~ 8 ). One possibility is to use Eq. (j25|) with m c equal to m T SSEH Better 
still, we numerically integrated the exact expressions for AC(/) provided in Ref. 1551 
for arbitrary values of my, obtaining a smooth dependence on the value of m c . For 
completeness we repeated this detailed analysis for all light fermions. As expected, 
the result depends very mildly on the values chosen for their masses. Employing 
the values m u = md — 0.3 GeV, m s = 0.5 GeV, m c = 1.5 GeV and m;, = 4.5 GeV, and 
adding to Eq. (|24|) the contribution of the remaining fermionic two-loop diagrams 
studied in Ref. 45, for M H = 150 GeV we obtain a™ {2 loop fer) = -4.68 x 10~ 8 . In 
this evaluation we also included the tiny O(10 -9 ) contribution of the j-Z mixing 
diagrams, suppressed by (1 — 4 sin 2 6* w ) ~ 0.1 for quarks and (1 — 4sin 2 6> w ) 2 for 
leptons, via the explicit formulae of Ref. [38] 

The sum of the fermionic and bosonic two-loop EW contributions described 
above gives a!p w (2 loop) = —7.74 x 10 -8 , a 14% reduction of the one-loop 
result. The leading- logarithm three- loop EW contributions to the muon q— 2 
were determined to be extremely small via renormalization-group analyses EHESI 
We assigned to our t lepton g— 2 EW result an additional uncertainty of 
0[a™(2 loop)(a/7r)ln(M 2 /m 2 )] - O(10" 9 ) to account for these neglected three- 
loop effects. Adding a™ (2 loop) to the one-loop value of Eq. (f2"3")l we get our total 
EW correction (for M^ISO GeV) 

a™ = 47.4(5) x 10~ 8 . (26) 

The uncertainty allows M H to range from 114 GeV up to ~300 GeV, and reflects 
the estimated errors induced by hadronic loop effects (m u and aid can vary between 
70 MeV and 400 MeV), neglected two- loop bosonic terms, and the missing three- 
loop contribution. It also includes the small errors due to the uncertainties in M top 
and m T . T he value in Eq. (f26|) is in agreement with the prediction a^ w = 47(1) x 
10 _8 JS3l25lwith a reduced uncertainty. As we mentioned in Sec. l3.lt the EW estimate 
of Ref. 33, af = 55.60(2) x 10~ 8 , mainly differs from Eq. |26]) in that it doesn't 
include the two-loop corrections. 

4. The Hadronic Contribution 

In this section we will analyze a^ AD , the contribution to the r anomalous magnetic 
moment arising from QED diagrams involving hadrons. Hadronic effects in (two- 
loop) EW contributions are already included in a™ (see the previous section). 
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4.1. Leading-order Hadronic Contribution 

Similarly to the case of the muon g— 2, the leading-order hadronic contribut ion to 
the t lepton anomalous magnetic moment is given by the dispersion integral!^ 

" ^ <x(°>(e+e- -^hadrons) K T (s) 



12tT 3 

where the kernel K T (s) is a bounded function of energy monotonously increasing to 
unity at s — > oo, and cr(°)(e + e~ — ► hadrons) is the total hadronic cross section of 
the e + e~ annihilation in the Born approximation. In Fig. [4] we plot the ratio of the 
kernels in the r lepton and muon case. Clearly, although the role of the low energies 
is still very important, the different structure of K T compared to K^, induced by 
the higher mass of the t, results in a relatively higher role of the larger energies. 




Fig. 4. Ratio of the kernels K T (s)/K fl (s). 



The history of these calculations is not as rich as that of the muon. The first 
calculation performed in 1978 in Ref. [55] was based on experimental data available 
at that time below 7.4 GeV, whereas at higher energies the asymptotic QCD pre- 
diction was used. Ten years later, a rough estimate was made in Ref. 1591 based on 
low energy e + e~ data. In Ref. [33]the contribution of the p meson was estimated by 
integrating the approximation obtained using the Breit-Wigner curve, while other 
contributions used the data. The accuracy of the calculation was considerably im- 
proved in Refs. [501 [52 where, below 40 GeV, only data were used. In Ref. [55J data 
were only used below 3 GeV (together with the experimental parameters of the 
J ftp and T family states). In our opinion this can significantly underestimate the 
resulting uncertainty. In addition, in the same reference, data from r lepton decays 
were extensively used; as it is known today, this leads to higher spectral functions 
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these calculations are summarized in Table [TJ 

Table 1. Calculations of the leading-order hadronic part of a T . 

Author Year a" LO x 10 s 

S. NarisonEHl 1978 370 ± 40 

B.C. Barish and R. Stroynowski ESI 1988 ~ 350 

M.A. Samuel et al. S3 1991 360 ± 32 

S. Eidelman and F. Jegerlehnerl 60 ! 61 ! 1995 338.4 ± 2.0 ±9.1 

S. NarisonESl 2001 353.6 ± 4.0 

This work 2007 337.5 ± 3.7 

M. Benmerrouche et al.lMI 1993 197-246 

F. Hamzeh and N.F. NasrallahED 1993 280 ± 20 



B. Holdom et al.ESl 1994 320 ± 10 

A.E. DorokhovE3 2005 310 ± 20 



For completeness, in the second part of Table [T] we also show purely theoretical 
estimates. The analysis based on QCD sum rules performed in Ref. [53] gives results 
which strongly depend on the choice of quark and gluon condensates. QCD sum rules 
are also used in Ref. [55] In Ref. [55J the authors use a nonlocal constituent quark 
model for the description of the photon vacuum polarization function ilhad(<7 2 ) at 
space-like momenta and obtain <z" LO = 3.2(1) x 10~ 6 , close to the estimates based 
on the experimental data. They also show that a simpler model with constituent 
quark masses independent of momentum is strongly dependent on the values chosen 
for the quark masses. For example, with m u — md = 330 MeV and m s = 550 MeV 
their result is a™° = 2.2(1) x 10~ 6 , i.e., significantly smaller than the previous 
estimate. They could reproduce the value a^ LO = 3.2(1) x 10 -6 using m„ = mj, = 
m s = 201 MeV. In a recent analysis using the instanton liquid model the author 
obtains (3.1 ±0.2) x 10 _6 P^ All these estimates somewhat undervalue the hadronic 
contribution and have rather large uncertainties. 

We updated the calculation of the leading-order contribution using the whole 
bulk of experimental data below 12 GeV, which include old d ata co mpiled in 
Refs.[5Dl[51[5Jl as well as the recent datasets from the CMD-J^MH] anc i SNe P° I 71 I 72 J 
experi ment s in Novosibirsk, a nd f rom the radiative return studies at KLOE in 
Frascat^Sl and BaBar at SLACKS] The improvement is particularly visible in the 
channel e + e _ — > 7r + 7r~, where four ne w i ndependent measurements exist in the 
most important p meson region: CMD-2,ES1 SNDpB and KLOrJ^Sl ( se e Fig. [5]). Our 
result is 

a* LO = 337.5(3.7) x 10~ 8 (28) 

(we recently presented a preliminary estimate of this value in Ref. I75|) . The break- 
down of the contributions of different energy regions as well as their relative fractions 
in the total leading-order contribution are given in Table [2] The contribution of the 
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Fig. 5. CMD-2, SND and KLOE data on |i<V| 2 in the p meson energy range. 

p meson energy range is still important, but its relative wei ght is s maller than in the 
case of the muon anomaly, 51.3% compared to about 72%!^21S31 The contributions 
of the narrow resonances (J/tp and T families) are included in the corresponding 
energy regions. It is worth noting that uncertainties of the contributions from the 
hadronic continuum are larger than that of the very precise 2?t one. The over all un- 
certainty is 2.5 times smaller than that of the previous data-based prediction! ^^ 



4.2. Higher-order Hadronic Contributions 

The hadronic higher-order (a 3 ) contribution a" HO can be divided into two parts: 
a™° = a™°(vp) + a™°(lbl). The first one is the 0(a 3 ) contribution of diagrams 
containing hadronic self-energy insertions in the photon propagators. It was deter- 
mined by Krause in 1996!^ 

< HO (vp) = 7.6(2) x 1(T 8 . (29) 

Note that naively rescaling the muon result by the factor / (as it was done 
in Ref. [33]) leads to the totally incorrect estimate a" HO (vp) = (-101 x 1CT 11 ) x 
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Table 2. Summary of the contribu- 
tions to a^ LO from various energy in- 
tervals. 



v/a, GeV 


Aa™° x 10 8 


Aafo, % 


2tt, < 2 


173.3 ± 1.6 


51.3 


u> 


15.0 ± 0.4 


4.4 


<t> 


21.1 ± 0.4 


6.3 


0.6-2.0 


48.2 ± 2.4 


14.3 


2.0-5.0 


57.9 ± 2.0 


17.2 


5.0-12.0 


16.9 ± 1.1 


5.0 


> 12.0 


5.1 


1.5 


Total 


337.5 ± 3.7 


100.0 



m^/mf l = —29 x 10 -8 (the a™°(vp) value is from Ref.[76j); even the sign is wrong! 

The second term, also of 0(a 3 ), is the hadronic light-by-light contribution. Sim- 
ilarly to the case of the muon g— 2, this term cannot be directly determined via a 
dispersion relation approach using data (unlike the leading-order hadronic contri- 
bution), and its evaluation therefore relies on specific models of low-energy hadronic 
interactions with elect romagnet ic currents. Actually, very few estimates of a™°(lbl) 
exist in the literatureJ^EHZSI anc [ a \\ f them were obtained simply rescaling the 
muon results a™°(lbl) by a factor m^/rrvt. Following this very naive procedure, 
the < HO (lbl) estimate varies between < HO (lbl) = [80(40) x lO" 11 ] x (m^/rn^) = 
23(11) x 10~ 8 , and < HO (lbl) = [136(25) x 10~ n ] x (m^/m*) = 38(7) x 10~ 8 , ac- 
cording to the values chosen for a™°(lbl) from Refs. 1771 and [751 respectively. 

These very naive estimates fall short of what is needed. Consider the function 
A^\mi/mj, lbl), the three-loop QED contribution to the g— 2 of a lepton of mass mi 
due to light-by- light diagrams involving loops of a fermion of mass m,j (see Sec. l2.2p . 
The exact expression of this function, computed in Ref. [30] for arbitrary values of 
the mass ratio mi /mj , is rather complicated, but series expansions were provided in 
the same article for the cases of physical relevance. In particular, if mj 3> mi, then 
A% (mi/mj, lbl) ~ (mi/mj) 2 . This implies that, for example, the (negligible) part 
of a™°(lbl) due to diagrams with a top-quark loop can be reasonably estimated 
simply rescaling the corresponding part of a" HO (lbl) by a factor m^/mi. On the 
other hand, to compute the dominant contributions to a" HO (lbl), i.e. those induced 
by the light quarks, we need the opposite case: mj <C mi = m T . In this limit, 
(mi /m 3 lbl) does not scale as (mi/mj) 2 , and a naive rescaling of a™°(lbl) by 
m^/m^ to derive a" HO (lbl) leads to an incorrect estimate. 

We therefore decided to perform a parton- level estimate of a" HO (lbl) based on 
the exact expression for A% {mi /mj , lbl) using the quark masses recently proposed 
in Ref. [7J5] for the determination of a™°(lbl): m u = m d = 176 MeV, m s = 305 
MeV, m c — 1.18 GeV and m& = 4 GeV (note that with these values the authors 
of Ref. [75] obtain a" HO (lbl) = 136 x 10~ n , in perfect agreement with the value in 
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Ref . 1751 - see also Ref . [5(71 for a similar earlier determination). We obtain 

a™°(lbl) = 5(3) x 1(T 8 . (30) 

This value is much lower than those obtained by simple rescaling of a™°(lbl) 
by m^/m^. The up-quark provides the dominating contribution; the uncertainty 
&™°(lbl) = 3 x 1CT 8 allows m u to range from 70 MeV up to 400 MeV. Further 
independent studies (following the approach of Ref. 78, for example) would provide 
an important check of this result. 

The total hadronic contribution to the anomalous magnetic moment of the r 
lepton can be immediately derived adding the values in Eqs. (|28l) . (f29f and ([30j) . 

fl HAD = a H LO + a HHo (vp) + qi^ ho (lbl) = 350.1 (4.8) x 10~ 8 . (31) 
Errors were added in quadrature. 



5. The Standard Model prediction for a T 

We can now add up all the contributions discussed in the previous sections to derive 
the SM prediction for a T : 

4 M = a? ED + a E T w + < LO + a™° (vp) + a™° (lbl) , (32) 

where 

a q ED = n7 324(2) xl0~ 8 

a™ =47.4(5) xl0~ 8 

< LO = 337.5 (3.7) xl0~ 8 

a™°(vp) =7.6(2) xlO- 8 

a™°(lbl) =5(3) xlO" 8 

(the sum of the hadronic contributions is given in Eq. (|3ip ). Adding errors in quadra- 
ture, our final result is 

a s T M = 117 721 (5) x lO" 8 . (33) 

The present PDG limit on the anomalous magnetic moment of the r lepton was 
derived in 2004 by the DELPHI collaboration from e + e~ — > e + e _ r + r _ total cross 
section measurements at ^fs between 183 and 208 GeV at LEP2P 

- 0.052 < a T < 0.013 (34) 

at 95% confidence level. The authors of Ref. [5] also quote their result in the form of 
central value and error: 

a T = -0.018(17). (35) 

Comparing this result with Eq. (|33jl (their difference is roughly one standard devi- 
ation), it is clear that the sensitivity of the best existing measurements is still more 
than an order of magnitude worse than needed. A reanalysis of various measure- 
ments of the cross section of the process e + e~ — > t + t~ , the transverse r polarization 
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and asymmetry at LEP and SLD, as well as of the decay width T(W — > tv t ) at LEP 
and Tevatron, allowed to set a stronger model-independent limitP^ 

- 0.007 < a T < 0.005. (36) 

Other limits on a T can be found in Refs. 82. 



6. Conclusions 

In this article we reviewed and updated the SM prediction of the r lepton g— 2. Up- 
dated QED and electroweak contributions were presented, together with new values 
of the leading-order hadronic term, based on the recent low energy e + e~ data from 
BaBar, CMD-2, KLOE and SND, and of the hadronic light-by-light contribution. 
These results were confronted in Sec. [5]to the available experimental bounds on the 
t lepton anomaly. 

As we already mentioned in the Introduction, quite generally, NP associated 
with a scale A is expected to modify the SM prediction of the anomalous magnetic 
moment of a lepton I of mass mi by a contribution af r ~ mf /A 2 . Therefore, given 
the large factor m^/m 2 ~ 283, the g— 2 of the r lepton is much more sensitive than 
the muon one to EW and NP loop effects that give contributions ~ mf, making its 
measurement an excellent opportunity to unveil (or just constrain) NP effects. 

Another interesting feature can be observed comparing the magnitude of the EW 
and hadronic contributions to the muon and r lepton g— 2. The EW contribution to 
the g—2 of the r is only a factor of seven smaller than the hadronic one, compared 
to a factor of 45 for the g—2 of the muon. Also, while the EW contribution to 
a^ M is only a factor of three larger than the present uncertainty of the hadronic 
contribution, this factor raises to 10 for the r lepton. If a NP contribution were 
of the same order of magnitude as the EW one, from a purely theoretical point 
of view, the g—2 of the r would provide a much cleaner test of the presence (or 
absence) of such NP effects than the muon one. Indeed, if this were the case, such 
a NP contribution to the r lepton g—2 would be much larger than the hadronic 
uncertainty, which is currently the limiting factor of the SM prediction. 

Unfortunately, the very short lifetime of the r lepton makes it very difficult to 
determine its anomalous magnetic moment by measuring its spin precession in the 
magnetic field, like in the muon g—2 experiment. Instead, experiments focused on 
high-precision measurements of the r lepton pair production in various high- ener gy 
processes, comparing the measured cross sections with the QED predictions! *^^ As 
we can see from Eq. ([34j> . the sensitivity of the best existing measurements is still 
more than an order of magnitude worse than that required to determine a T . 

Nonetheless, the possibility to improve such a measurement is certainly not 
excluded. For example, it was suggested to determine the r lepton g factor taking 
advantage of the radiation amplitude zero whic h o ccurs at the high-energy end 
of the lepton distribution in radiative r decays!^ This method requires a very 
good energy resolution and could perhaps be employed at a r-charm or B factory 
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also benefiting from the possibility to collect very high statistics. It is not clear 
whether the huge data samples at B factories will result in a corresponding gain 
for the limits on a T . Indeed, LEP measurements were rather limited by systematic 
uncertainties, which were of the order of 2-3% for the discussed processes and, until 
now, experiments at B factories have not yet reached such a level of accuracy in 
the absolute measurements of the tot al cr oss sections. However, a search for the r 
lepton electric dipole moment at Belle^l showed that with the appropriate choice 
of observables, using full information about events, the improvement in sensitivity 
can be proportional to the square root of luminosity, i.e., determined mainly by 
statistics. One can hope that this is also the case with the determination of a T . 
A similar method to study a T using radiative W decays and potentially very high 
data samples at LHC was suggested in Ref. |85] Yet another method would use 
the channeling in a bent crystal similarly to the suggestion for the measurement 
of magnetic moments of short-living baryons!^ This method has been successfully 
tested by the E761 collaboration at Fermilab, which measured the magnetic moment 
of the S + hyperon!S3 In the case of the r lepton, it was suggested to use the 
decay B + — > t + v t , which would produce polarized r leptons!^ In 1991, when 
this suggestion was published, the idea seemed completely unlikely. However, in the 
era of B factories, when the decay B + — > t + v t is already observed by the Belle 
collaboration,^ anc i the possibility of a Super- B factory is actively discussed, this 
is no longer a dream. Even more promising could be the realization of this idea 
in a dedicated experiment at a hadron collider with its huge number of B mesons 
produced and a more suitable geometry. We believe that a detailed feasibility study 
of such an experiment, as well as further attempts to improve the accuracy of the 
theoretical prediction for a T , are quite timely. 
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